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Abstract- A third order parabolic operator £g, typical of a non linear wave 
operator Cq perturbed by viscous terms, is analyzed. Some particular solutions 
J^ ' related to £o a-i'e explicitly determined and the initial value problem for /C^ is 

^Sl ' considered. The parabolic-hyperbolic behaviour is analyzed and rigorous approx- 

t^^ , imations of the solution are achieved. 

cn 

^^ ' Riassunto- Si prende in esame un operatore parabolico del terzo ordine C^ 

tipico di perturbazioni introdotte su equazioni iperboliche non lineari. Si ricavano 
anzitutto classi di soluzioni dell'equazione non pertubata e si considera per C^ il 
problema di valori iniziali in tutto lo spazio. Mediante una opportuna equazione 
integrale si determinano stime utili alia valutazione rigorosa di approssimazioni di 



C^ . strato limite. 
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Introduction 



In a vast number of realistic mathematical models (ecology, supercon- 
ductivity, spread of epidemics, neurobiology, viscoelasticity) , the evolution 
is characterized by deep interactions between wave propagation and diffu- 
sion. These models are generally given by non linear hyperbolic equations 
perturbed by viscous terms described by higher - order derivatives with 
small diffusion coefficients e. Such small parameters imply very long times 
for the transmission of signals, at least when the diffusion is the principle 
process involved. 

In order to have an enough knowledge of the dynamic response, the 
order of diffusion effects and the time - intervals where the wave behavior 
prevails must be estimated. The valuation of these aspects leads to the 
analysis of non linear parabolic - hyperbolic boundary layers. A typical 
example is related to equations such as: 

(1.1) Ce u = [edxxt + c^ dxx - dtt - adt]u = f{x,t,u,Ux,ut) 

where u = u{x,t) and e, a, c are positive constants. 

The parabolic equation (1.1) includes models of the viscoelasticity [l]-[3] 
and biology [4]. In particular, when / = 7 +senu, one has the perturbed 
sine - Gordon equation (PSGE) which is basic in superconductivity [5], [6]. 

In this paper, the initial value problem Ve for the PSGE in all of the 
space is considered and the fundamental solution Kg{x,t) of the linear 
operator Ce is explicitly determined. So, when / is linear (/ = f{x,t)), 
the problem Vg is quite solved. More, K,, is a C°^ function endowed with 
properties like those of the fundamental solution of the heat operator [7]. 
For this, when / is not linear, it is convenient to reduce the problem Vg to 
an integral equation whose kernel is just K,, . 

The qualitative analysis of this equation allows not only existence and 
uniqueness results, but rigorous estimates when e — )• too. Let Cq the 
operator to which Cs is reduced for e = and let w the solution of the 
following initial value problem Vq : 



(1.2) Co w = f{x,t,u) with w{x,0) = fo{x), wtix,0) = /i(x). 
In the superconductive case, it results [8] 

(1.3) f{u) — f{w) = senw — senu 

and the error v = u — w related to the approximation of u by means 
of the wave solution w, is the solution of a problem Ve with initial data 
vanishing and a source term given by 

(1.4) Fw{v) = senw — senu — ewxxt- 

Then, according to several classes of special initial data (/o, /i,),the 
solution w of the reduced problem Vq is explicitly determined and, in 
corrispondence, the remainder term v = u — w is rigorously estimated, 
together with the time - interval where the diffusion effects are neglegible 
for e — 7- 0. As meamingful result one can see that whatever a positive 
constant A; < 1 is prefixed, the evolution is characterized by the travelling 
wave w and by diffusion effects of the order of e , in each interval- time 
[0, Te], with Ts defined in (6.12). 



2 Statement of the problem 

If T is a positive constant and 

nT = {{x,t) -.xeR, 0<t<T}, 
let consider the following initial value problem Ve : 



dxx{£Ut + c^u) - dtiut + au) = f{x,t,u), ix,t) £ Qt, 



(2.1) 



u{x,0) = fo{x), utix,0) = fi{x), 



X £ R, 



where /, /o, /i are arbitrary specified functions, and e, a, c, are positive 
constants. 

When e = 0, the parabolic equation (2.1)i turns into the hyperbohc 
telegraph equation 



(2.2) 



(c^ dxx - dtt - adt)w = f{x, t, w) 



and the problem Vg changes into a problem Vq for ■w{x, t) which has the 
same initial conditions of Vg. 

When e is neglegible, a possible boundary - layer region could appear for 

T = cx). 

In order to estimate the influence of the dissipative term e Uxxt on the 
wave behavior of w , the difference 



(2.3) 



u{x, t, e) — w{x, t) = v{x, t, e) 



is to be evaluated. For this, one has the following problem Vv related to 
the remainder term v : 



£Vxxt + c^Vxx-vtt-avt = Fw{x,t,v) (x,t)G (7t, 
(2-^) \ v{x,0)=0, vt{x,0)=0, xeR, 



where the source term F^ is given by: 



(2.5) 



F^{x,t,v) = f{x,t,w + v) - f{x,t,w) - ewxxt- 



3 Linear case : fundamental solution and proper- 
ties 



When Fyu = F{x,t) is linear and the Laplace transform is applied to 
the problem (2.4), the transform v{x,s) of the solution v{x,t) is: 

(3.1) Hx, s)= f k,{x - i, s) Fit s) dt 

JR 

where F is the transform of F and 



- \x\ ^s(s+a)/(es+c2) 

(3-2) K,{x,s)= —====^=. 

To obtain the inverse Laplace tranform C~[ of the function (3.2), let 
r = \x\/e, b = (? je and 

^-r^s(s+a)/(s+6) 

(3.3) G{t,s) 



2\/^V(s + «)(s + 6)' 
so that: 



(3.4) k, (x, s)= ^ G(r, >,)=£,[ ^ * G(r, t) 

If I„(-z) denotes the modified Bessel function, let 



t -^ 
(3.5) G(r,t) = -^ / ^ e-^^*-'^) /o(rV'(& "«)(*- «)/t' d^- 



4 a/tt^ Jo ^ v^ 



Then, the following theorem holds: 

Theorem 3.1- For all r > 0, the Laplace integral Ct G{r,t) converges 
absolutely in the half-plane Kes > max{—a, —b), and one has: 

(3.6) /:tG{r,t) = G{r,s) , 

with G and G defined by (3.3)-(3.5). 

Proof - By Fubini- Tonelli theorem it results 



{3.7)CtG=-^ e'^'-+^^^ e-^'+'^'Io{rJ{b-a)z/v)dz 

4:^/-Ke Jo Vyjv Jq "^ 



and further [9] it is: 



(3.8) r e-(*+^)^ /o ( r^{b-a)z/v ) dz 



^ s + b 

By means [9] of known formulae, (3.7) and (3.8) imply (3.6). The absolute 
convergence holds whatever the constants a, b may be. In fact, when 
a > b, for all real z > 0, one has: 



(3.9) \Mr^{a-b)z\ < 1 < Io{r^{b - a) z ). 



Remark 3.1 - The function G(r,t) given by (3.5) is defined also when 
r = \x\ / ^/e = 0. In fact, if in the integral of (3.5) one puts : r^ {t—v)/v = 
z, one has 



a 



(3.10) G(M)= j^ r .-*-^" 'A^^ ,, 

4 V Tret JO V-z + ^ 

and for r = 0, it results: 



e'^* f^ ^-^ ^ , r, — . dz 



(3.11) G(0,t) = -^= / e-t /o(^(6-a)z) 



4V7ret Jo ^ \/^ 

1 °+f'/ ^ , b — a , 



As consequence, by (3.4) it follows 



(3.12) K,{x,t) = f G{r,T) 

Jo 



dT 



/o v^7r(t - r) 

and, for x = , it is: 



(3.13) i^,(0,t) =1 f f " ^ loi^r) dr. 

2 Jo y^TTe{t-T) 2 



Moreover, it's easy to prove the following theorem: 

Theorem 3.2- The function K^ defined by (3.10)- (3.12) is a C°°{9.t) 

solution of the equation Lg^u = 0. More, when a <b = (? /e, K^ is never 
negative in VLt- I 

As (3.10)- (3.12) show, the fundamental solution K^{x, ■ ), qua function of 
X, is an even function which results positive when a < b and so 

/•oo roo 

(3.14) 0< / Ke{x-(,t)d( =2 Ke{y,t)dy. 

Jo Jo 

More by (3.2) one has 



oo . 1 I 



(3.15) 0<2/ K,{y,s)dy = ^ = £,[_(i_e 

JO s"^ + as a 

and the absolute convergence of the integral (3.6) implies the following 
estimate: 



Theorem 3.3- When a < b it results 



rco 1 

(3.16) 0< / Ke{x-i,t)di = {-){!- e~''*) <l/a. 

Jo a 



This property of K,; is useful to various estimates for the non linear 
case. 



4 Reduced equation and travelling waves 



As for the model of Superconductivity, the approximation w is charac- 
terized by the source / = senw -|- 7, so that one has 



(4.1) Wxx - wtt - awt = smw +7 

where it's assumed (? = 1 without loss of generality. 

Let ijj an arbitrary function and let n(V') the Lobachevsky's angle of 
parallelism [10], defined by 

{n(V') = 2 arctg e"^ = 2arctg e'"^ [V' > 0] 
n(v) = vr - n(-v) [^<o], 

and such that 

(4.3) senli{il)) 



cosh{ip) 



A class of travelling wave solutions of (4.1) can be explicitly evaluated 
whatever 7 may be; it suffices to consider functions ip = V'(0 such as: 



(4.4) w = U[ij(0] = ^(x,t), with S, = -(x-t). 

a 

Then (4.1) becomes 

(4.5) — a(pt = simp +7 
and the function ip is given by 

1 + ^ cosnyip) 

According to all the values of 7 , indicating by k an arbitrary constant 
of integration, one has the following special solutions of the non linear 
equation (4.1). 

- Case 7 = 

(4.7) w{x,t) = n(^ + k), 

- Case 7 = 1. If ^ — A; = ^, one has: 



(4.8) w{x,t) = 2arctg{-^^] 



Case 7^ < 1. If a = \J\ — 7^, {a / k) ^ = ^, one has: 



, a ,1 + k e^ 1 
(4.9) w{x,t) = 2 arctg[ - {- 



7 1 - k e^ a 



-Case 72 > 1. If CT = v/7^ - 1, {a/2)S.=S., one has 



(4.10) w{x,t) =2arctg{ - [l + atg{S, + k)]}. 

7 



5 An integral equation 

Consider now the non linear problem defined in (2.4) and assume that 
the prefixed approximation w{x, t) is compatible with the following assump- 
tions on the function F^ given by (2.5): 

Conditions A 

1) Fw{x,t,v) is defined and continuous on the set 

T = {{x,t,v) G r^T X (-00, +oo)} 

2) For each A; > and for \v\ < k , the function Fu,{x,t,v) is uni- 
formly Holder continuous in x and t for each compact subset of 
Qr • 

3) There exists a constant Cp such that 

I F^ (x, t, vi) - Fy, (x, t, V2) \ < Cf\vi-V2\ 
holds for all 1^1,1^2 and w. 

4) Fw is bounded in 17^ for bounded v and w. 

Definition 5.1 -For each prefixed w{x,t), by a solution of the prob- 
lem Vv we mean a function v continuous and bounded in Qt such that, 
whatever e > may be, the derivatives vt, vu, Vxx, Vxxt are continuous in 
r^T and verify (2.4). 

According to the formal representation (3.1) it results: 

(5.1) V{x,t)= f dT f Ke{\x-C\,t-T) F^[C,T,v{C,T)]dC 

Jo Jr 
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and the conditions A imply that, for a solution f (x, t) of the problem 
"P^jthe function X{S,,t) = Fw[6., t, i^{S,,t)] is bounded in Qt and is uni- 
formly Holder continuous in S^ and continuous in r on each compact of 

If 



(5.2) \\Fw\\t = sup sup |A(,^,t)|, 

o<T<t x<aR 



by (5.1) and theorem 3.3 one deduces 



(5.3) \vix,t)\ < (t/a) II F^ lit xeR, t>0 

and so, for each compact subset of — oo < x < oo, < t < T, the 
function v is uniformly Holder continuous. Consequently, by Conditions 
A, the function F^ [ x, t, v{x, t) ] is uniformly Holder continuous too. More, 
by means of the properties of the kernel K^, it's possible the following 
statement. 

Theorem 5.1 - For each prefixed w, the initial-value problem (2.4) 
admits a unique solution if and only if the integral equation (5.1) possesses 
a unique solution v continuous and bounded for {x,t) G (— oo,+oo) x 

(o,r). I 

Consider now a positive rj < T and let 



(5.4) \\v\\r] = sup sup I v{x,t) 

0<t<»7 xeR 



(5.5) Bfj = {v{x,t) : V £ C{ {—oo, cx)) x [0, 77] ) and \\v ||^ < 00 } 
The set Br/ is a Banach space and the mapping 



(5.6) J='v{x,t) 



f dr f K,{x-tt-T)F^[^,T,v{tT)]dt 
Jo Jr 
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owing to theorems 3.2 - 3.3, maps Bri into Brj. Then, by means of the well- 
known techniques of the fixed point, the following theorem can be proved. 

Theorem 5.2 - When Fu,{x,t,v) satisfies the conditions A, then the 
initial- value problem Vv given by (2.4) admits a unique solution. I 

In conclusion, according to the prefixed initial data fo{x), fi{x) of the 
problem Ve, one can solve explictly the reduced problem Vq by means 
of the research of travelling waves outlined in sect. 4 and by numerous class 
of solutions of sine- Gordon equations available in literature. 

When the solution w of Vq is determined, the source term F^ is 
known and its properties can be analyzed in order to show that theorems 
5.1 and 5.2 can be applied. Then, in this case, the integral equation (5.1) 
allows to estimate the remainder term v. 



6 Estimate of the diffusion - An example 



For example, consider the approximation w induced by the initial con- 
ditions 



(6.1) 



w{x,0) = 2arctg{e^^"-) 



and let £, = - {x — t). Then, the problem Vq given by (4.1) - (6.1)admits 
the solution 



(6.2) w{x,t) = 2 arctg{e^) 

which represents a travelling wave with speed c = 1 . It results: 



2 e"^^ - 6 e^^ + 1 

(6.3) w^^^ = - - e^ ———— (^ei?; 
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and so 

(6.4) I w^^t I < /3 V (x, t) G J^T 

where the constant /3 depends only on a . As consequence, the source 
F,uj (x, t, v) related to the problem Vv for the superconductive model is 

(6.5) Fw{x,t,v) = sen {v + w) — senw — ewxxt 
and satisfies the conditions A. More, by (6.4)one has: 

(6.6) \ F^{x,t,v)\ < \v{x,t,e)\ +ef3. 

Then, by theorems 5.1 and 5.2, one deduces that the error v related to 
the approximation Us ^ w is such that: 



-^.7)\v{x,t,e)\ < f dT f [\v{tT,e)\ + e ^ ]\ K,{x - ^, t - T)\dC 

Jo JR 



2 

When e — 7- 0, it results a < b = — and Kg; > according to theorem 
3.2. If one puts 



(6.8) reit) = sup \v{x,t,e)\, 

R 

by (6.7) and theorem 3.3 it follows 

(6.9) < re{t) < - / re{T)dT + - et 

a Jo a 

and so the Gronwall Lemma implies 

(6.10) < re{t) < [/3(r/a) e^/"] e VtG[0,r]. 

13 



This estimate allows to specify the infinite time - intervals where the 
effects of diffusion are of the order e^ with k < 1. In fact, whatever 
< k < 1 may be, let 

Then, for T < T^, by (6.10) it results 

(6.12) < rs{t) < /3 e'^^^" e < e^ Vt G [O.Te]. 



So, when e is vanishing, the evolution of the superconductive model is 
characterized by the travelling wave w given in (6.2) and by diffusion effects 
which are of the order of e^ in each interval- time [0, T^], with T^ given 
by (6.11) and k < 1. 
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